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ABSTRACT 

We analyze the 2-point correlation function (2PCF) of galaxy groups identified from 
the 2-degree Field Galaxy Redshift Survey with the halo based group finder recently 
developed by Yang et al. (2004b). With this group catalogue we are able to estimate 
the 2PCFs for systems ranging from isolated galaxies to rich clusters of galaxies. The 
real-space correlation length obtained for these systems ranges from ~ 4 /i~^Mpc to 
~ 15/i~^Mpc, respectively. The observed correlation amplitude (and the corresponding 
bias factor) as a function of group abundance is well reproduced by associating galaxy 
groups with dark matter haloes in the standard ACDM model. Redshift distortions 
are clearly detected in the redshift-space correlation function, the degree of which is 
consistent with the assumption of gravitational clustering and halo bias in the cosmic 
density field. In agreement with previous studies we find a strong increase of the 
correlation length with the mean inter-group separation. Although well determined 
observationally, we show that current theoretical predictions are not yet accurate 
enough to allow for stringent constraints on cosmological parameters. Finally, we use 
our results to explore the power-law nature of the 2PCF of galaxies. We split the 
2PCF in 1-group and 2-group terms, equivalent to the 1-halo and 2-halo terms in halo 
occupation models, and show that the power-law form of the 2PCF is broken, when 
only including galaxies in the more massive systems. 
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1 INTRODUCTION 

In the standard cold dark matter (CDM) cosmogony galax- 
ies are assumed to form in virialized dark matter haloes. 
Theoretically, the properties of the halo population can be 
studied in great detail with the use of high-resolution A*'- 
body simulations and sophisticated analytical models. Ob- 
servationally, however, dark matter haloes can only be de- 
tected indirectly; either through their gravitational lensing 
of background sources, or by using galaxies and/or X-ray 
gas as tracers of the dark matter potential wells. In this pa- 
per we investigate the clustering of dark matter haloes using 
the second method. 

Based on galaxy kinematics. X-ray studies and gravita- 
tional lensing effects, it is now well established that clusters 
of galaxies are associated with the most massive dark mat- 
ter haloes. Observations show that clusters of galaxies are 
strongly clustered. The cluster-cluster two-point correlation 
function, £,cc{r), is roughly a power law, £,cc{r) = [ro/r)" , 
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with Q ~ 1.8 and with a correlation length, tq, that is much 
larger than that of galaxies (see Bahcall 1988 for a review; 
Croft et al. 1997; Park & Lee 1998; Bahcall et al. 2003). 
This is in good agreement with clusters being associated 
with massive dark matter haloes, which are expected to be 
strongly clustered from the fact that they are associated 
with high peaks in the initial density field (e.g.. Kaiser 1984). 

The correlation length of clusters is also observed to 



increase with the mean inter-cluster separation d 



-1/3 



with n the number density of objects (e.g. Bahcall & West 
1992). Since richer clusters (where richness expresses the 
number of galaxy members) are rarer objects, this relation 
between ro and d is equivalent to a relation between tq and 
cluster richness. Associating richer clusters with more mas- 
sive haloes, this, again, is in good agreement with theoret- 
ical predictions. Mo, Jing & White (1996), using the halo 
bias model developed by Mo & White (1996), found that 
the observed r^-d relation can be well described in terms 
of halo-halo correlation functions in the CDM cosmogony. 
In addition, these authors showed that the relation between 
correlation length and mean separation can be used to con- 
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strain models of structure formation. In particular, they 
showed that the observed ro-d relation is well reproduced 
by a ACDM model with n^fi = 0.3 and h = 0.7, but dif- 
fers significantly from the predictions for a CDM model with 
f^m,o = 1 and h = 0.5 (see also Bahcall et al. 2003). 

The ro-d relation for poorer systems can be probed by 
studying the correlation function of galaxy groups (Zandi- 
varez et al. 2003; Padilla et al. 2004). Using groups of galax- 
ies selected from the 2dFGRS (Eke et al. 2004), Padilla et 
al. (2004) found that the so-d relation (with so the correla- 
tion length in redshift space) obeyed by clusters extends to 
poor groups. Unfortunately, the connection between galaxy 
groups and dark matter haloes is less straightforward as for 
rich clusters, simply because the smaller number of galax- 
ies involved in individual groups makes it harder to identify 
systems that are physically associated. The correspondence 
between galaxy groups and dark matter haloes may there- 
fore depend significantly on the group finder used, compli- 
cating the interpretation of the observational results. In or- 
der to overcome this problem, one needs a group finder that 
associates galaxies according to their common dark matter 
haloes. 

In a recent paper (Yang et al. 2004b, hereafter YMBJ), 
we developed a halo-based group finder that is optimized 
for grouping galaxies that reside in the same dark matter 
halo. We tested the performance of this group finder exten- 
sively using mock galaxy redshift surveys constructed from 
the conditional luminosity function model (Yang, Mo, van 
den Bosch 2003; van den Bosch, Yang, Mo 2003; Yang et al. 
2004a), and found that our group finder is more successful 
than the conventional friends-of-friends (POP) algorithm in 
associating galaxies according to their common dark matter 
haloes. In particular, our group finder performs also reliably 
for very poor systems, including isolated galaxies in small 
mass haloes. 

In this paper we analyze the 2-point correlation func- 
tion (2PCP) of the galaxy groups identified by YMBJ. As 
we will show, our group catalogue allows us to determine 
2PCFs for vastly different systems, ranging from isolated 
galaxies to rich clusters. Using detailed mock galaxy redshift 
surveys (hereafter MGRSs), we show that the group corre- 
lation functions are closely related to those of dark matter 
haloes. The paper is organized as follows. In Section 2, we 
briefly describe our group finder, and summarize the proper- 
ties of the group catalogues obtained from the 2dPGRS and 
the MGRSs. In Section 3 we estimate the 2PCP of galaxy 
groups in the 2dPGRS. The relation between the correlation 
of galaxy groups and that of dark matter haloes is exam- 
ined in Section 4. In Section 5, we analyze how the corre- 
lation length of galaxy groups depends on the abundance 
of the systems in consideration, and we compare the obser- 
vational results with theoretical predictions in Section 6. In 
Section 7, we use our results to discuss how one can under- 
stand the galaxy-galaxy correlation function in terms of the 
group-group correlation function and the galaxy occupation 
in groups. Finally, in Section 8, we summarize our results. 

Unless stated otherwise, we consider a fiat ACDM 
cosmology with flm ~ 0.3, f^A ~ 0.7 and h — 
Ho/{WO kms"^ Mpc"^) = 0.7 and with initial density fluc- 
tuations described by a scale-invariant power spectrum with 
normalization as — 0.9. All distances are calculated using 
this cosmology. 



2 GALAXY GROUPS AND DARK MATTER 
HALOES 

2.1 The Group Finder 

In a recent study (YMBJ), we developed a halo-based group 
finder that can successfully assign galaxies into groups ac- 
cording to their common haloes. The basic idea behind our 
group finder is similar to that of the matched filter algo- 
rithm developed by Postman et al. (1996) (see also Kep- 
ner et al. 1999; White & Kochanek 2002; Kim et al. 2002; 
Kochanek et al. 2003; van den Bosch et al. 2004a,b), al- 
though we also made use of the galaxy kinematics. In sum- 
mary (see YMBJ for details), the group finder starts with 
an assumed mass-to-light ratio to assign a tentative mass to 
each potential group. This mass is used to estimate the size 
and velocity dispersion of the underlying halo that hosts the 
group, which in turn is used to determine group member- 
ship (in redshift space). This procedure is iterated until no 
further changes occur in group memberships (see Appendix 
for more details). We tested the performance of our group 
finder in terms of completeness of true members and contam- 
ination by interlopers, using detailed MGRSs. The average 
completeness of individual groups is ^ 90 percent and with 
only ~ 20 percent interlopers. Furthermore, the resulting 
group catalogue is insensitive to the initial assumption of 
the mass-to-light ratios, and the group finder is more suc- 
cessful than the conventional POP method (Eke et al. 2004) 
in associating galaxies according to their common dark mat- 
ter haloes. 



2.2 2dF Groups and Mock Catalogues 

In YMBJ we applied the group finder described above to 
the final public data release of the 2dPGRS. This observa- 
tional sample contains about 250, 000 galaxies with redshifts 
and is complete to an extinction-corrected apparent magni- 
tude of bj « 19.45 (CoUess et al. 2001). The survey volume 
of the 2dFGRS consists of two separate declination strips in 
the North Galactic Pole (NGP) and the South Galactic Pole 
(SGP), respectively, together with 100 2-degree fields spread 
randomly in the southern Galactic hemisphere. When identi- 
fying galaxy groups, we restricted ourselves only to galaxies 
with redshifts 0.01 < z< 0.20 in the NGP and SGP regions. 
Only galaxies with redshift quality parameter g > 3 and with 
redshift completeness > 0.8 were used. This left a grand to- 
tal of 151, 820 galaxies with a sky coverage of 1124 deg^. We 
obtained a group catalogue of 78, 708 systems which in total 
contain 104, 912 galaxies. Among these systems, 7251 are 
binaries, 2343 are triplets, and 2502 are systems with four 
or more members. The vast majority of the groups (66, 612 
systems) in our catalogue, however, consist of only a single 
member. Note that some faint galaxies are not assigned to 
any systems, because it is difficult to decide whether they 
are the satellite galaxies of larger systems, or whether they 
are the central galaxies of small haloes. 

As discussed in YMBJ, it is not reliable to estimate 
the (total) group luminosity based on the assumption that 
the galaxy luminosity function in groups is similar to that 
of field galaxies. We therefore used a more empirical ap- 
proach to estimate the group luminosity Lis, defined as 
the total luminosity of all group members brighter than 
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Figure 1. The rcdshift distributions of galaxy groups for tliree different bins in Lig (as indicated). Open squares with errorbars are tfie 
mean and 1-cr variance of the number counts for groups in 8 independent MGRSs, while solid dots correspond to the number counts of 
groups in the 2dFGRS. Solid lines indicate the expected value for a constant group number density. As shown in YMBJ, groups with 



Mb J — 51og/i = —18. In the Appendix we describe in de- 
tail how Li8 is estimated for each group. As demonstrated 
in detail in YMBJ, Lig is tightly correlated with the mass 
of the dark matter halo hosting the group, and can be used 
to rank galaxy groups according to halo masses. Fig 1 plots 
the redshift distributions of groups detected in the 2dFGRS 
(sohd dots) and in our MGRSs (open circles) . The solid line 
corresponds to a constant number density, and is shown for 
comparison. As already shown in YMBJ, the group cata- 
logue is virtually complete over the entire redshift range 
(0.01 <z< 0.20) for groups with Lxa > W^^'^h'^LQ (right- 
hand panel). For groups with smaller Lis, the catalogue is 
incomplete: groups with 10^°/i"^ L© < Lis < W^^'^h'^ L© 
are only complete to 2; ~ 0.13 (middle panel), while those 
with 10^-^ h^'^ Lq < ii8 < 10^° h^'^ Lq are complete down 
to z ~ 0.08 (left-hand panel). Note that the redshift dis- 
tributions of the MGRS agree nicely with the 2dFGRS, in- 
dicating that we have properly accounted for the various 
incompleteness effects when constructing our mock surveys 
(see Yang et al. 2004a and van den Bosch 2004a for details) . 
With these considerations, we can construct volume-limited 
group samples by ranking all groups according to their Lis. 
The brightest A'^ groups then form a volume- limited subsam- 
ple. Using this ranking-technique we construct 9 subsamples 
Oi, where i = 1, 2, • • •, 9 correspond to different choices of A^ 
and the maximum redshift Zmax- Rather than characterizing 
different subsamples by A'^ and Zmax , we use the mean group 
separation, d — n^^'^, where n is the number density of 
groups in the subsample. Table 1 lists the subsamples thus 
selected, and which form the observational data base for our 
analyses. 

In YMBJ, we also applied our group finder to eight 
MGRSs constructed using exactly the same selection criteria 
as the 2dFGRS (see Yang et al. 2004a for details). Note that 
in the present MGRSs, we also include fiber collisions in 
a way as described in van den Bosch et al. (2004b). Here 
we make use of these mock group catalogues to test the 
relation between the groups and the dark matter haloes. 
For this purpose, we generate eight dark halo catalogues 
from the eight MGRSs. In which, we select all dark matter 



haloes in our 'virtual universe' with 0.01 < z < 0.20 that 
are within the area of the sky covered by 2dFGRS where the 
completeness is larger than 0.8. Note that these haloes are 
not exactly the same as those corresponding to all selected 
groups, because the later are not complete due to the survey 
selection effect. Subsamples of mock groups and dark matter 
haloes are constructed in the same way as the 2dF group 
samples, i.e. according to the Lis ranking for mock groups, 
or according to halo mass ranking for dark matter haloes. 
We denote these two sets of subsamples as Mi (for mock 
groups) and Hi (for dark matter haloes). Subsamples Oi, 
Mi and Hi all have the same number of objects for a given 
i. Since we have 8 independent MGRSs, for each i we have 
8 independent mock subsamples and 8 halo subsamples. All 
errorbars quoted below are based on the scatter among these 
subsamples. 



3 THE TWO-POINT CORRELATION 
FUNCTION 

We compute the group-group (or halo-halo) 2PCF ^{rp,-K) 
using the following estimator 



ar 



(RR) (DP) 
{DRY 



(1) 



with (DD), (RR), and (DR) the number of group-group, 
random-random, and group-random pairs with separation 
(rp,7r) (Hamilton 1993). Here Vp and n are respectively the 
pair separations perpendicular and parallel to the line-of- 
sight. Explicitly, for a pair (81,82), with 8i = cZiii/Ha, we 
define 

TT = — , Tp = Vs ■ S - 7r2 (2) 

Here 1 = (l/2)(si + S2) is the line of sight intersecting the 
pair, and s = si — 82. 

Except for 01, all 2dF samples listed in Table 1 contain 
sufficient numbers of galaxy groups for a proper determina- 
tion of the 2PCF. Fig 2 shows the contour-plots for the 
^(rp, tt) of some of these samples. Panels from upper left to 
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Table 1. The 2dFGRS group correlation functions 



Sample 


N 


•2^max 


d 
Mpc/h 


soi-y = 1.8) 
Mpc/h 


ro(7 = 1.8) 
Mpc/h 


''0 

Mpc/h 


7 


b/b(06) 


/3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Ol 


250 


0.20 


43.77 


19.30 ±1.73 


- 


- 


- 


- 


- 


02 


500 


0.20 


34.74 


16.09 ±1.42 


15.05 ±0.83 


15.44 


2.38 


1.79 ±0.11 


0.22 ±0.27 


03 


1000 


0.20 


27.57 


14.84 ±0.63 


12.79 ±0.58 


13.53 


2.16 


1.58 ±0.07 


0.17 ±0.16 


04 


2000 


0.20 


21.88 


12.68 ± 0.41 


11.83 ±0.48 


11.82 


1.79 


1.34 ±0.04 


0.28 ±0.15 


05 


4000 


0.20 


17.37 


10.84 ± 0.32 


9.62 ±0.32 


9.36 


1.72 


1.19 ±0.03 


0.36 ±0.07 


06 


8000 


0.20 


13.79 


9.26 ±0.22 


8.11 ±0.17 


8.02 


1.77 


1.00 


0.35 ± 0.08 


07 


16000 


0.20 


10.94 


8.12 ±0.17 


6.94 ±0.23 


6.54 


1.68 


0.87 ±0.02 


0.48 ± 0.07 


08 


16000 


0.13 


7.22 


6.24 ±0.20 


4.77 ±0.12 


4.78 


1.86 


0.63 ±0.03 


0.50 ±0.12 


09 


8000 


0.08 


5.67 


4.55 ±0.40 


3.55 ±0.23 


3.70 


1.85 


0.48 ± 0.05 


0.68 ±0.26 



Column (1) indicates the sample ID. Columns (2) (number of groups) and (3) (redshift range: 0.01 < z < 



^) indicate the selection 



criteria. Column (4) lists the mean separation {d = n~^'^) of the selected groups. Columns (5) and (6) list the redshift space and 
real-space correlation lengths, respectively, obtained fitting t;{s) with ^(s) = (s/so)^'®(5 < s < 15 h~^Mpc) and Wp{rp) using ^(r) = 
(r/ro)"'^'*(3 < rp < 15 h~^Mpc), Column (7) and (8) indicate the ro and 7 obtained fitting Wp{rp) using g(r) = (r/ro)"'. Column (9) 
indicates the bias of groups relative to that of the fiducial 06 sample, and column (10), finally, lists the 13 parameter. All the errbars 
listed in this table are 1-a variances obtained from the scatter among 8 mock group samples. 
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Figure 2. The two-point correlation function, ^(fp, tt), for various group samples (as indicated) extracted from the 2dFGRS. From upper 
left to upper right, then lower left to lower right, the samples have smaller mean inter-group separations d (see Table 1), indicating an 
increased inclusion of less massive systems. Note that samples with smaller d reveal a more pronounced flattening of the contours (see 
also Section 5 and the right-hand panel of Fig. 7). 
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Figure 3. The projected correlation function, Wp{rp) (upper panels), and redsliift space correlation function, ^(s) (lower panels), of 
groups and dark matter haloes. Solid dots correspond to the groups extracted from the 2dFGRS, while solid and dashed lines indicate 
the same results but obtained for the samples of mock galaxies and their corresponding dark matter haloes, respectively. The errorbars 
associated with the solid lines indicate the I-ct variance obtained from the eight independent MGRSs. Results are shown for three different 
values of the mean group separation, d, as indicated. 



upper right and from lower left to lower right correspond to 
samples 03 - 08. Notice that these £,{rp,-K) look very dif- 
ferent from those of galaxies (e.g., Hawkins et al. 2003): the 
only deviation from isotropy is a flattening of the contours 
at large separations due to the infall motion induced by the 
gravitational action of large scale structure. Unlike for galax- 
ies, no finger-of-God effect on small scales is present, due to 
the fact that groups themselves are virialized objects rather 
than test particles in larger virialized potentials. As we will 
see in Section 5, this absence of virial motions on small scales 
makes the interpretation of the redshift distortion easier. 

Since the redshift-space distortion only affects tt, the 
projection of ^(rp,7r) along the n axis can get rid of the 
infall induced distortions and give a function that is more 
closely related to the real-space correlation function. This 
projected 2PCF, 'Wp{rp), is related to the real-space 2PCF, 
£,{r), through a simple Abel transform 



Wp(rp)= / ^(rp,7r)d7r = 2 / ^(r) 



r dr 



yjr'^ - rl 



(3) 



(Davis & Peebles 1983). Therefore, if the real-space 2PCF 
is a power-law, ^(r) = (ro/r)^, the projected 2PCF can be 
written as 



£(7/2-1/2) /n, 



(4) 



The black dots in the upper panels of Fig 3 show the 
projected correlation function Wp{rp) of 2dFGRS groups 
estimated from £,{rp,ii) using eq. (3) with the integration 
range set to |7r| < 40 h,^^Mpc. For comparison, we also plot 
Wp{rp) for the mock groups (solid line with errorbars). The 
three panels correspond to samples with mean separation 
d = 34.74 /i~^Mpc, 17.37 fe^^Mpc, and 7.22 /i"^Mpc, as in- 
dicated. Overall, the agreement between data and mock is 
extremely good. An exception is the large scales in sam- 
ple 08 (d = 7.22 /i~^Mpc)T, where the Wp{rp) of the mock 
groups is significantly underestimated. This is due to the fact 
that this sample occupies a small, nearby volume, which in 
our MGRSs is represented by a small box-size simulation 
that does not properly sample the large(r) scale structure 
(see Yang et al. 2004a for details). 

In most previous studies of group-group correlation 
functions, the redshift-space 2PCF ^(s), rather than the 
real-space 2PCF, was used to represent the clustering 
strength (Croft et al. 1997; Park & Lee 1998; Zandivarez 
et al. 2003; BahcaU et al. 2003; Padilla et al. 2004), In order 
to allow for a comparison we also compute the redshift-space 
correlation functions which are shown in the lower panels of 



T We find a similar discrepancy between data and model for sam- 
ple 09, which is limited to an even smaller volume, with an even 
smaller mean separation, than sample 08, 
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Fig. 3. Here again, the results for the mock samples match 
those of the 2dFGRS samples remarkably well, except at 
large radii in sample 08. 



4 THE RELATION BETWEEN GALAXY 
GROUPS AND DARK MATTER HALOES 

So far we have focused on the 2PCFs for groups in the 2dF- 
GRS and in our MGRSs. We now examine whether these 
results can be understood in terms of 2PCFs between dark 
matter haloes in the ACDM concordance cosmology. Since 
the clustering properties of CDM haloes are well understood 
(Mo & White 1996; 2002; Sheth & Tormen 1999; Jing 1998; 
Sheth, Mo & Tormen 2001; Jenkins et al. 2001; Seljak & 
Warren 2004), such a connection between the populations 
of galaxy groups and dark matter haloes enables us to un- 
derstand the clustering of groups in a cosmological context. 

As mentioned above, the luminosity of a group, Lis, is 
tightly correlated to the mass of its host halo. Therefore, 
groups ranked by the value of Lis may be used to represent 
dark matter haloes ranked by halo mass. To check this, we 
compare the correlation functions of mock group samples 
(Ml, M2, etc) with those of dark matter halo samples (HI, 
H2, etc). The results are shown as solid and dashed lines, 
respectively, in Fig. 3. Note that the correlation function of 
mass-ranked dark matter haloes matches that of Lig-ranked 
groups remarkably well. 

In order to facilitate a more qualitative comparison, we 
fit Wp{rp) with a single power-law of the form (4) over the 
range 3 h~^Mpc < rp < 15 h~^Mpc. The goodness of fit 
is based on a simple x^ criterion, where the errors used for 
each data point are obtained from the scatter among 8 in- 
dependent mock samples (the errors due to cosmic variance 
are typically larger than the statistical errors on each indi- 
vidual Wp{rp) measurement). Over the Vp range considered 
here, a power law is an acceptable model. We treat the slope 
7 either as a free parameter or keep it fixed at a value of 
7 = 1.8. In the latter case, the fit is used to determine only 
the correlation length ro. The difi'erences in the correlation 
lengths estimated with fixed or free 7 is less than 10%. Fig. 4 
shows ro (obtained keeping 7 fixed) as a function of mean 
group separation for both mock groups and dark matter 
haloes. The agreement between the groups and dark mat- 
ter haloes is remarkably good, especially for massive/bright 
systems (note that they have been offset from each other 
by Alogd = 0.03 for clarity). At small values of d, i.e. for 
faint groups and low-mass haloes, the groups have slightly 
lower correlation lengths than the dark matter haloes. This 
discrepancy is at least partly due to the incompleteness of 
the 2dFGRS (which we have mimicked in our MGRSs) . Due 
to this incompleteness, which is not present for the dark 
matter haloes, the true mean inter-group separation is over- 
estimated. This effect is less important for larger groups; 
although some of the member galaxies are missed, they still 
contain sufficient members to be identified as a group. 

We have also estimated the redshift-space correlation 
lengths. So. In this case, we adopt a simple power- law model, 
5(s) = (s/so)^**, to fit ^(s) over the range 5 h~^Mpc < s < 
15 h^^Mpc. The lower limit of s adopted here is larger than 
that of rp used in fitting Wp{rp), because the redshift-space 
correlation is not well described by a power law at smaller 




Mpc) 



Figure 4. Relation between correlation length, ro, and mean 
inter-group separation, d, for groups (solid circles) and dark mat- 
ter haloes (solid squares) in the MGRSs (errorbars indicate the 
1-(T scatter among the eight independent mock catalogues). For 
clarity, the results for the dark matter haloes have been shifted 
to the right by Alogd = 0.03. Note the good agreement between 
groups and haloes, indicating that groups ranked by luminosity 
can be compared directly to dark matter haloes ranked by halo 
mass. Thick and thin solid lines correspond to theoretical pre- 
dictions based on the halo bias models of SW04 and SMTOl, 
respectively. Note that the difference between these two model 
predictions is larger than the scatter among our eight MGRSs 
(sec Section 6 for a detailed discussion). 



separations (see lower panels of Fig. 3). As for the real- 
space correlation lengths, we find extremely good agreement 
between the so of mock groups and dark matter haloes (not 
shown) . 

All these results provide strong support for a tied cor- 
relation between group luminosity and halo mass, clearly 
demonstrating that the groups ranked by luminosity can 
be compared meaningfully to dark matter haloes ranked by 
halo mass. 



5 ABUNDANCE DEPENDENCE OF GROUP 
CORRELATION FUNCTION 

Having established a tied correlation between group lumi- 
nosity and halo mass, we now return to our 2dFGRS group 
catalogue. Using the same fitting procedure as described 
above, we determine the correlation lengths ro and Sq as 
well as the slope 7 for 2dFGRS groups. Results are listed in 
Table 1 (columns 5 to 8). In Fig. 5 we plot the correlation 
lengths ro (left panel) and So (right panel) as a function of 
d, obtained assuming a fixed slope of 7 = 1.8. For compari- 
son, the right-hand panel also shows the SDSS results (open 
squares) obtained by Bahcall et al. (2003) and the 2PIGG 
results (open circles) obtained by Padilla et al. (2004) from 
the 2dFGRS. All three measurements are in excellent agree- 
ment with each other. Fitting our ro-d and so-d relations by 
power laws, we obtain ro = 1.11 d°'^^ and so = 1.88 d°'^^. 
These power laws are shown as dot-dashed lines in Fig 5. 
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Figure 5. Left panel: The relation between rg and d for groups selected from the 2dFGRS (solid dots). Errorbars indicate the 1-cr variance 
from 8 independent mock group samples. The solid and dashed lines are model predictions for a ACDM cosmology with erg = 0.9 and 
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and 2PIGG (Padilla et al. 2004). The dot-dashed line corresponds to the best-fit power-law, so = 1.88 dO^^. 



The correlation strength of a sample can also be 
described by the ratio of its projected correlation func- 
tion and that of a fiducial sample. Fig 6 plots the ratio 
Wp{rp)/wp^06{'''p) as function of Vp, where 06 has been used 
as the fiducial sample. Note that Wp(rp)/wp,o6(rp) is roughly 
constant with Vp, indicating that the slope of the correlation 
function is roughly the same for different samples. The over- 
all amplitude of the ratio, however, has a clear trend with 
d. To illustrate this we define the relative bias 6/6(06) as 
the mean value of the ratio Wp{rp)/Wpfie{i'p) in the range 
5 /i"^Mpc < rp < 10 h'^Mpc, and plot b/b{06) as function 
of mean group separation in the left-hand panel of Fig. 7 
(numerical values are listed in column (9) of Table 1). 

For a given sample the value of so is systematically 
larger than ro (cf. columns 5 and 6). This is due to the 
enhancement of clustering in redshift-space due to gravita- 
tional infall. To quantify this redshift distortion, we use the 
model of Kaiser (1987; see also Hamilton 1992). According 
to linear theory, the infall velocities around density pertur- 
bations affect the observed correlation function as 

Ciin(rp,7r) = Co(s)Po(m) +C2(s)P2(m) +C4(*-)^4(m) , (5) 

where Vi{^) is the /*'' Legendre polynomial, and ^ is the 
cosine of the angle between the line-of-sight and the redshift- 
space separation s. According to linear perturbation theory 
the angular moments can be written as 



Co(s) 



i + ^ + ^lCM 



us) = {^ + ^)[ar)~m] 



with 



ar) + l^r) - liir) 



(6) 

(7) 
(8) 



^W = ^y C(r>"dr', (9) 

and 

iir) = -^ f\{r')r"dr' . (10) 

^ Jo 

In the above expressions, /3 is the linear distortion param- 
eter, which can be written as /3 = O^^/fe, where b is the 
bias parameter of the objects under consideration. Given the 
real-space correlation function ^(r) — (r/ro)'', which can be 
obtained from Wp{rp), eq. (5) can be used to model (,{rp,Tv) 
on linear scales. By comparing the model predictions with 
the observed ^(rp,7r), one can easily obtain the value of /3. 
We use a simple x^ fit of the observed C('"p) ^) in the range 
8fe~^Mpc < s < 20/i~^Mpc to eq. (5) to probe the only free 
parameter l3. In the fitting, each data point for ^(rp,7r) is 
weighted by the error based on the scatter among 8 indepen- 
dent mock samples. The right-hand panel of Fig. 7 plots the 
P values thus obtained for both mock groups (squares with 
errorbars) and 2dFGRS groups (solid dots) . Although there 
is significant scatter, there is a clear trend that groups with a 
smaller mean separation d (i.e., less luminous groups) have 
larger /3 and thus a stronger distorted redshift-space cor- 
relation function (this is also directly visible from Fig. 2). 
The numerical values of /3 for the 2dFGRS groups are listed 
in column (10) of Table 1, together with the 1-a variances 
obtained from the scatter among the 8 MGRSs. 



6 COMPARISON WITH THEORETICAL 
PREDICTIONS 

The tests described in the previous sections show that the 
abundance-dependence of the group-group correlation func- 
tion can be explained in terms of the halo-halo correlation 
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function (after all, we constructed our group finder to asso- 
ciate galaxies according to tiieir common darlc matter iialo) . 
This suggests that we may compare the 2dFGRS group- 
group 2PCF with halo-halo correlation functions predicted 
by current models of structure formation in order to con- 
strain cosmological parameters. 

The mean number density of dark matter haloes with 
mass M > Mi can be estimated through. 



n{M)dM, 



(11) 



where n{M) is the mass function of dark matter haloes, 
which can be estimated analytically from the Press- 
Schechter formalism (e.g. Press & Schechter 1974; Sheth, 
Mo & Tormen 2001, hereafter SMTOl). The mean bias for 
haloes with mass exceeding Mi can be estimated from 



b{M > Ml) = - 



n{M > Ml 



n{M)b{M)dM , 



(12) 



Ml 



where 6(M) is the bias parameter of dark matter haloes (Mo 
& White 1996; Jing 1998; Sheth & Tormen 1999; SMTOl; 
Seljak & Warren 2004, hereafter SW04). Throughout we use 
the halo mass function of SMTOl, which has been shown to 
be in excellent agreement with numerical simulations (e.g., 
Jenkins et al. 2001; White 2002). For the halo bias parame- 
ter, we use the models of both SMTOl and SW04 for com- 
parison. 

Using d = n~^''^ and ro — b^^^'^ ro.DM, with ro.DM the 
correlation length of the dark matter, we compute ro{d). 
Here the linear power spectrum is computed using the trans- 
fer function of Eisenstein & Hu (1998), which properly ac- 
counts for the baryons, while the non-linear power spectrum, 
which is required in calculating the dark matter correlation 
function and ro,DM, is computed using the fitting formula 
of Smith et al. (2003). The sohd lines in Fig 4 show the 
model predictions thus obtained using the bias models of 
SW04 (thick hue) and of SMTOl (thin line). The difference 
between the two bias models is quite large, and much larger 
than the l-o variance among our 8 MGRSs. Thus, although 
the ro-d relation can now be accurately determined from ob- 
servational data the current models for halo bias are not yet 
accurate enough to allow one to obtain stringent constraints 
on model parameters (see below). 
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In the left panel of Fig 5 we compare the results ob- 
tained for the 2dFGRS groups with various theoretical pre- 
dictions. The thin solid line corresponds to a standard 
ACDM model with ag = 0.9, obtained using the bias model 
of SMTOl. The thin dashed line indicates the prediction 
for the same bias model but with erg = 0.7. Based on this, 
one might conclude that the observational data is in bet- 
ter agreement with erg = 0.9. However, if we use the bias 
model of SW04 (thick solid and dashed lines), the ag = 0.7 
cosmology matches the data better. Note that although the 
bias model of SW04 may be more accurate than earlier mod- 
els, the uncertainty of the bias parameter at the massive end 
is still 10 to 20 per cent, which is much larger than the error 
on the observational results. Clearly, the halo bias model 
has to be improved further, in order to make full use of the 
constraining power of the present observational results. 

Fig. 7 compares the theoretical predictions for the rel- 
ative bias &/6(06) (left panel) and the redshift-distortion 
parameter /? = O!^ /h (right panel, assuming Q.m ~ 0.3) 
with our observational results from the 2dFGRS. Due to the 
normalization at a given d, the predicted relation between 
the relative bias and d is quite similar for different models 
(i.e., the relative bias is fairly insensitive to the value of ug). 
More importantly, all model predictions are in good agree- 
ment with the observational data. The observed value of j3 
as a function of d is also well described by the theory, but 
the errorbars are too big to provide stringent constraint on 
model parameters. 



7 UNDERSTANDING THE SHAPE OF THE 
CORRELATION FUNCTION OF GALAXIES 

It is well known that the real-space correlation function of 
(normal) galaxies is remarkably well described by a single 
power law for r ;$ 10 /I'^Mpc. Given that the mass corre- 
lation function predicted for typical ACDM cosmologies is 
significantly curved on these scales, it is important to under- 
stand the origin of this power-law behavior. Jing et al. (1998) 
were the first to show that , if the number of galaxies in a dark 
matter halo increases with halo mass as a power law, with 
a power index moderately below unity, and if the number 
density distribution of galaxies in massive haloes has approx- 
imately the same profile as the dark matter, the observed 
power-law shape of the galaxy correlation function can be 
reproduced. This kind of galaxy bias on small scales is now 
well understood in the current halo occupation model (Pea- 
cock & Smith 2000; Seljak 2000; Scoccimarro et al. , 2001; 
Jing, Borner & Suto 2002; Berlind & Weinberg 2002; Bul- 
lock, Wechsler & Somerville 2002; Scranton 2002; Berlind et 
al. , 2003; Yang, Mo & van den Bosch 2003; van den Bosch, 
Yang & Mo 2003). In the halo model, the 2PCF of galaxies 
can be decomposed into two terms: 



C(r)=Cih(r)+6h(r), 



(13) 



where ^ih represents the correlation due to pairs of galaxies 
within the same halo (the "1-halo" term) , and ^211 describes 
the correlation due to galaxies that occupy different haloes 
(the "2- halo" term). In the standard ACDM model ^(r) has 
a characteristic scale at r ~ 1 - 2 /i~^Mpc, where the domi- 
nating contribution to the 2PCF makes a transition from the 
1-halo term to the 2-halo term. Therefore, some departure 



from a pure power law is expected for populations of galaxies 
for which the 1-halo and 2-halo terms are not well balanced. 
In fact, such a departure, albeit small, has recently been 
found in the projected correlation function of SDSS galaxies 
(Zehavi et al. 2004). 

The various analyses in the previous sections have 
shown that our groups selected from the 2dFGRS are nicely 
related to dark matter haloes. Therefore, we can directly 
measure the '1-halo' and '2-halo' terms of the 2PCF, by 
simply determining whether both galaxies of a pair reside 
in the same group (this pair adds to the 1-group term), or 
whether they reside in two different groups (in which case 
the pair adds to the 2-group term). Fig. 8 plots the pro- 
jected two-point correlation functions of 2dFGRS galaxies 
that are associated with groups of different abundances. In 
addition to Wp,tot(''p), we also indicate the '1-halo' (1-group) 
and '2-halo' (2-group) terms as dot-dashed and dashed lines, 
respectively. Each galaxy pair is weighted by Wg = l/(ciCj), 
where Ci is the survey completeness at the position of galaxy 
'i'. On scales Vp ^ 3/i^^Mpc the projected correlation func- 
tion mp,tot(fp) is dominated by the '2-halo' term, while on 
smaller scales (r^ ^ 1/i^^Mpc) the '1-halo' term dominates. 
Note also that for galaxies residing in massive systems (i.e., 
large d), the projected correlation function clearly deviates 
from a pure power-law; for these galaxies, the '1-halo' term is 
significantly enhanced. This occurs because the '1-halo' term 
contribution from a large group is proportional to Ng{Ng — l) 
(with Ng the number of galaxies in the group), while the 
'2-halo' term contribution is proportional to b^, which in- 
creases with halo mass at a slower rate. When adding more 
galaxies to the sample hosted by smaller haloes (i.e., de- 
creasing d), Wp^tot{rp) becomes better described by a pure 
power-law. In particular, as shown by the solid line in the 
lower right-hand panel, when all galaxies (151,820 galaxies 
with 0.01 < 2 < 0.20 and completeness > 0.8) are included, 
the correlation function is well represented by a single power 
law. 

As an additional test, we estimate Wp,tot(rp) for a sam- 
ple in which we remove all galaxies in most luminous groups 
such that the total number of galaxies in the sample is 
halved. Thus, this sample contains only galaxies in low-mass 
haloes. The correlation function for this sample is shown as 
the long dashed line in the lower right panel of Fig 8. As for 
the complete sample, the projected correlation function of 
this sample is well described by a power law, with a slightly 
shallower slope as for the complete sample. 

Therefore, we conclude that the 2PCF of 2dF galaxies 
reveals a power-law form as long as sufficiently many small 
mass groups (halos) are included. For galaxies hosted by 
massive haloes, however, the 2PCF can deviate significantly 
from a pure power law. 



8 CONCLUSIONS 

We have measured the 2PCFs for galaxy groups in the 2dF- 
GRS group catalogue constructed by YMBJ using a halo- 
based group finder. We have shown that the current data al- 
lows one to estimate the correlation function accurately for a 
wide range of different systems, ranging from isolated galax- 
ies to rich clusters of galaxies. Ranking groups according to 
their luminosities, Lig, we have studied how the correlation 
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Figure 8. The projected two-point correlation function Wp{rp) of galaxies in different samples of groups. The solid lines indicate the 
total correlation functions, Wp^tot{rp), of all galaxies in group samples 04 (upper left), 06 (upper right), 08 (lower left) and of all 
(151,820) galaxies in the 2dFGRS sample with 0.01 < z < 0.20 and completeness > 0.8 (lower right). The dot-dashed and short dashed 
lines indicate the corresponding '1-halo' terms, Wpn^(rp), and '2-halo' terms, uip_2h(^p)) respectively. Finally, the long dashed line in the 
lower right panel indicates the projected correlation function of that half of all galaxies that is not associated with the luminous groups. 
See text for details and discussion. 



of groups depends on group abundance. Consistent with pre- 
vious studies (e.g. Bahcall et al. 2003; Padilla et al. 2004), 
we found that the ampUtude of the correlation function in- 
creases with group luminosity (richness) . The dependence of 
the redshift-space correlation length so on the mean inter- 
group separation d can be quantified as so = l.SSdP'^^, 
while the real-space correlation length ro reveals a some- 
what steeper dependence: ro = 1.11 d"'^^. 

Using mock group catalogues, obtained from detailed 
mock galaxy redshift surveys, and the corresponding cata- 
logues of dark matter haloes, we have shown that the cor- 
relation functions of the 2dFGRS groups can be understood 
in terms of halo-halo clustering. The observed correlation 
length (and the corresponding bias factor) as a function of 
group abundance is well reproduced by associating galaxy 
groups with dark matter haloes in the standard ACDM 
model. In particular, the groups ranked by Lig match ex- 
tremely well with dark matter haloes ranked by mass. We 
found, however, that current theoretical predictions for the 
halo-halo correlation functions are not yet accurate enough 



to allow us to use the observational results to put stringent 
constraints on model parameters in the ACDM cosmogony. 
Analyzing the correlation function for galaxies associ- 
ated with different groups, we were able to bisect the 2PCF 
of galaxies in terms of a group-group correlation function 
and a term due to galaxies in the same group. Since our 
groups are closely related to dark matter haloes, this split 
correponds to the 1-halo and 2-halo terms of the correlation 
function. We have shown how the power-law form of the 
(projected) correlation function is broken when only consid- 
ering galaxies in massive halos, and how the balance between 
the '1-halo' and '2-halo' terms changes with halo mass. 
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APPENDIX A: THE HALO-BASED GROUP 
FINDER AND THE ASSIGNMENT OF Lis 

Al The group finder 

In a recent paper, Yang et al. (2004b) developed a halo- 
based group finder that can successfully assign galaxies into 
groups according to their common haloes. For completeness, 
we present a brief description of the group finder here, but 
refer the reader to Yang et al. (2004b) for details. 

The halo-based group finder consists of the following 
main steps: 

Step 1: Two different methods are combined to iden- 
tify the centres (and members) of potential groups. First we 
use the traditional Friends-Of-Friends (FOF) algorithm with 
very small linking lengths to assign galaxies into groups. The 
geometrical centres of all FOF groups thus identified with 
more than 2 galaxies are considered as centres of potential 
groups. Next, from all galaxies not yet linked together by 
these FOF groups, we select bright, relatively isolated galax- 
ies which we also associate with the centres (and members) 
of potential groups. 

Step 2: We estimate the luminosity of a selected po- 
tential group using 



U 



^ fciU) 



(Al) 



where Li is the luminosity of each galaxy in the group, and 
fc is the incompleteness of the survey. The total luminosity 
of the group is approximated by 



Lt. 



= L^ 



/p°° L(^{L)dL 
' J^ L(f,[L)dL 



(A2) 



where Lum is the minimum luminosity of a galaxy that can 
be observed at the redshift of the group, and (j){L) is the 
galaxy luminosity function. 

Step 3: From Ltotai and a model for the group mass-to- 
light ratio, wc compute an estimate of the halo mass associ- 
ated with the group in consideration. From this estimate we 
also compute the halo radius rigo, the virial radius rvir, and 
the virial velocity V^,\i — (GM/rvir)^ . The line-of-sight ve- 
locity dispersion of the galaxies within the dark matter halo 
is assumed to be a = Vvir/\/2. 

Step 4: Once we have a group centre, and a tentative 
estimate of the group size, mass, and velocity dispersion, we 
can assign galaxies to this group according to the properties 
of the associated halos. If we assume that the phase-space 
distribution of galaxies follows that of the dark matter parti- 
cles, the number density contrast of galaxies in redshift space 
around the group centre {= centre of dark matter halo) at 



redshift 2:group can be written as 



PAf(i?,A2) 



Ho S(i?) 



p(A2) 



(A3) 



Here Az = z — Zgroup and S(i?) is the projected surface den- 
sity of a (spherical) NFW halo, while the function p[Az)AAz 
describes the redshift distribution of galaxies within the 
halo. See Yang et al. (2004b) for the functional forms of 
E(7?) and p{/^z) used. 

Thus defined, Pm{R,^z) is the three-dimensional den- 
sity contrast in redshift space. In order to decide whether a 
galaxy should be assigned to a particular group we proceed 



as follows. For each galaxy we loop over all groups, and com- 
pute the corresponding distance {R, Az) between galaxy and 
group centre. Here R is the projected distance at the redshift 
of the group. If Pm {R, Az) > B, with _B = 10 an appropri- 
ately chosen background level, the galaxy is assigned to the 
group. If a galaxy can be assigned to more than one group, 
it is only assigned to the group for which Pm {R, Az) has the 
highest value. Finally, if all members of two groups can be 
assigned to one group according to the above criterion, the 
two groups are merged into a single group. 

Step 5: Using the group members thus selected we re- 
compute the group-centre and go back to Step 2, iterat- 
ing until there is no further change in the memberships of 
groups. Note that, unlike the traditional FOF method, this 
group finder also identifies groups with only one member. 



A2 The assignment of Lig to Groups 

As discussed in Yang et al. (2004b), it is not reliable to 
estimate the (total) group luminosity based on the assump- 
tion that the galaxy luminosity function in groups is similar 
to that of field galaxies. Therefore we used a more empiri- 
cal approach to estimate the group luminosity Lis, defined 
as the total luminosity of all group members brighter than 
Mtj - 5 log /i 
as follows. 



-18. The assignment of Lis to a group goes 



• We estimate the group luminosity Lgroup with eq. (Al) 
using only galaxies with Mbj — 5\ogh < —18.0. 

• We compute the absolute magnitude limit Miyj^um — 
5 log h at the redshift of the group in consideration. 

51og/i > —18.0, then we set L\s = 



If Mbj^Um 

,up; otherwise L\s 



,Mb 



igroup X f(Lgroup,Mtj,iim), wherc 
is the correction factor determined from 



groups at lower redshifts where the galaxy sample is com- 
plete down to Mb J - 5 log ft < -18.0. 

To determine the correction factor between Lgroup and 
Lis we first select all groups with z < 0.09, which cor- 
responds to the redshift for which a galaxy with Mbj — 
51og/i — —18 has an apparent magnitude equal to the 
mean limiting magnitude of the 2dFGRS (bj < 19.3). 
By applying further absolute magnitude limit cuts, we es- 
timate I/group (Afft J > Mb J, cut), and consider the ratio 
Lis/Lgroup as a function of Mbj,cut) (see Fig. 9 in Yang et 
al. 2004b). We fit the Lis / Lg^oup-Mbj ^cut relation to a func- 
tional form [log(Lgroup/ao)l"''^^'^'+"' '^*^, where AM = 
Mbj,iim — 5 log ft — 18.0, by adjusting ao, ai and a2, and 
use this fit result as an estimate for the correction factor, 
/(igroup, Mbj^iim)- Notc that the correction factor obtained 
in this way is an average, and is not expected to be accurate 
for individual groups. However, as demonstrated in detail in 
Yang et al. (2004b), Lis so defined is quite tightly correlated 
with the mass of the dark matter halo hosting the group in 
the mock catalogue, and can be used to rank galaxy groups 
according to halo masses. 
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